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Abstract 

A novel imaging principle based on the interaction of electromagnetic waves with a 
beam of relativistic electrons is proposed. Wave-particle interaction is assumed to 
take place in a small spatial domain, so that each electron is only briefly accelerated 
by the incident field. In the one-dimensional case the spatial distribution of the 
source density can be directly observed in the temporal spectrum of the scattered 
field. Whereas, in the two-dimensional case the relation between the source and the 
spectrum is shown to be approximately the Radon transform. 
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1 Introduction 



Telescopes, antenna arrays, radars, and imaging interferometers share a com- 
mon resolution criterion: the spatial extent of their aperture. Higher resolution 
means a larger lens, a larger reflector, or a longer interferometer base. At the 
present stage any significant improvement of the resolution is associated with 
extreme technical difficulties and costs, and for some applications we are ap- 
proaching a limit in this respect. Hence, it makes sense to consider alternative 
imaging principles, especially ones where the attainable resolution would not 
be entirely controlled by the spatial extent of some 'aperture'. 
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In a variety of cases both the observer and the observed object are in rel- 
ative motion with respect to each other. A well-known phenomenon related 
to this motion is the Doppler effect. It is widely exploited in radar [1], and 
acoustical imaging of blood vessels [2] , [3] , and is sometimes used to improve 
the resolution in astronomy [4] . Not to mention the famous red shift of stars 
and galaxies around us, which is the main source of cosmological information. 
However, in all of these imaging techniques the observer is actually considered 
to be at rest. Often, this is just a question of reference frame, and we simply 
find it convenient to relate ourselves with the one at rest. If for some reason 
the motion of the observer cannot be neglected in this way, then our natural 
intention is to compensate for it [5], [6]. 

In this paper we investigate the possibility of exploiting the observer's motion 
in electromagnetic imaging, rather than neglecting or correcting for it. Ap- 
parently, for this purpose we need an observer that moves with a relativistic 
velocity. This, however, presents not a problem, since an electron, which is 
easily accelerated to such velocities, is also the most natural "observer" for 
the electromagnetic field. An obvious, but not at all unique way to organize 
an imaging system of this kind is to let a beam of relativistic electrons inter- 
act with the fields radiated by some remote spatially inhomogeneous source 
distribution, say a group of stars. Contrary to standard (stationary) imaging 
instruments, here it is the observer's velocity, not the spatial extent of the 
aperture, that determines the attainable resolution. 

One of the problems in [7] (p. 249), where a plane monochromatic wave in- 
teracts with a single relativistic electron, contains a hint, which has inspired 
the present paper. In [7] it is claimed that upon scattering another plane wave 
can be observed in the far-field, which has a different frequency with respect 
to the incident one, depending on the angles of incidence and observation. 
Here we simulate an imaging setup and, therefore, consider incident fields 
due to a remote but spatially extended source. These fields interact with a 
short segment of a relativistic current. The short-segment assumption allows 
to circumvent the difficult problem of charge dynamics. Straightforward cal- 
culations presented below show that the relation between the spatial density 
distribution of the source and the (temporal) spectrum of the secondary field 
radiated by the electrons is approximately a Radon transform. This transform 
is well known in imaging theory and is amenable to invertion [8] . Mathemati- 
caly, the particular form of the Radon transform obtained here is very similar 
to the Doppler and Doppler-Radon transforms, which are extensively studied 
in the (non- relativistic) acoustical imaging of fiuid fiows [2] , [3] . 

The author does not claim to have covered all physical and mathematical 
aspects of this interesting problem. In fact, to arrive at the Radon transform 
in its simplest form we make several approximations, the significance of which 
must be further investigated both theoretically and experimentally. 



2 



2 Theory 



Let there be a remote spatiaUy inhomogeneous and time varying source of the 
electromagnetic field. We consider the scattering of this field by a segment 
of relativistic electric current, e.g. a beam of electrons moving with a rela- 
tivistic velocity in a particle accelerator or even a simple cathode ray tube. 
Neglecting the mutual interaction between the electrons, quantum effects, as 
well as the radiation reaction force, we would like to calculate the secondary 
electromagnetic field, scattered by this relativistic current into the far-field 
zone. 

Let Xe(t) denote the location of a single moving electron and let Dg be a 
bounded region of space occupied by an external source of the electromagnetic 
field (see Fig. 1). Presume that Dg is far enough both from the origin of 
the coordinate system and location Xe(i), so that the distance factor can be 
approximated by 

|Xs - Xe(te)| ~ |Xs| - ■ Xe(ie), for Xg G D^, (l) 

where ng = Xs/|xs|. In Gaussian units the incident electromagnetic field at the 
location of an electron is given by 

E-(x.,g = l / n.xK><y(^.'.)) dv;. (2) 

C J Xg 

H-(x„g = l / °.xy(^.*^) dv;. (3) 

C J Xg 

where J(x, t) is the electric current density in the remote source. Time ts in 
the right-hand-sides of (2)-(3) is retarded with respect to time te- 

The acceleration of a single electron due to these fields is given by the well- 
known relativistic formula [7] 



+/3nvXff°(xe,g" 



(4) 



where (3 = |ve|/c, ny = Ve/|ve|, and Vg is the electron velocity. This formula 
clearly shows that we need to take into account the action of both the electric 
and the magnetic fields of the source. Subsequently, an electron radiates into 
the far- field the following electric field [7] : 
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Fig. 1. Imaging setup with a relativistic detector (electron beam). 

Here, Xq denotes the location of the observation point. Time in the right- 
hand-side of this expression is retarded with respect to time to- Upon substi- 
tuting (2)-(3) in (4), and the result in (5), we arrive at 



E(Xo,to) 



c^me|xo| [1 - /? (no • ny)f 

lis X (lis X dt,3{y^,ts)) 

|Xs| 



[no - /?nv] X 



dK- 



(3'^ny n^ • J 

J xJ 



(6) 



dK 



dK 
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This expression describes the scattering by a single electron. The electric field 
scattered by a beam of non-interacting electrons is simply the sum of the 
fields scattered by each electron individually. All these fields can be computed 
using (6), provided that the electron velocities are given. In principle, these 
velocities are solutions of the relativistic equation of motion (4), which is a 
nonlinear equation and therefore very difficTilt. if not impossible, to solve. For 
our purposes, however, it is not necessary to know the velocities exactly. We 
rather need a reliable estimate on their variation with respect to the initial 
relativistic velocity. 

To obtain such an estimate we use the Taylor expansion 

Ve{to + At) ^ Ve{to) + dtVe\t=to^t + higher order terms, (7) 

and observe that Ve(t) can be approximated by the initial velocity up to the 
order At - the time of the field-particle interaction. Note, if At is sufficiently 
small, then the initial velocity approximation is suitable even for nonvanishing 
accelerations, and is, therefore, consistent with the idea of the secondary (i.e. 
acceleration related) radiation. Since At is approximately the time of flight 
of electrons across the interaction area, one can decrease At by decreasing 
this area. This can be achieved by either actually creating a relatively short 
relativistic current (radio frequencies) or by focusing the incident fields (op- 
tics) , so that the latter are different from zero only within a certain finite and 
relatively small spatial domain. Under these assumtions we may take in (6) 

Ve(te) ~ Ve(to) — COUStaut, 

f3{te) ~ f3{to) — constant, (8) 
ny{te) ~ nv(to) = constant. 

As far as the electron location Xe(te) is concerned, we are only interested in 
Xe(te) within the (small) area of interaction. For each electron this location 
depends on time linearly up to 0{At). Under these assumptions the distance 
expression (1) becomes 

|Xs| - lis • Xe(te) ~ |Xs| - te (Hs " Ve) + Ilg • Xe(to), (9) 

where Xe(to) is the initial location of the electron inside the area. 

We set the origin of the coordinate system to be within the area of interaction. 
Recall that the times at the source and the electron locations are retarded with 
respect to each other and with respect to the time at the observation point. 
In the initial velocity approximation these times are 

4 = ie - '"^^^^^ ~ + /3(ns • n.)] - ^, (10) 

c c 

t^^t,- l^'o-Mte)] ^^^^ ^^^^^^ . _ W _ ^11) 
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up to the order At. The arbitrary initial locations of electrons vanish, and, 
therefore, approximately the same time is now associated with all electrons 
inside the interaction area. From (lO)-(ll) we obtain the following relation 
between the times at the source location and the observation point: 

V c y 1 - /3(no • nv) c 



A particularly simple expression in the case of almost constant velocities is 
obtained for the following choice of the observation direction: 



n„ = 1. 



(13) 



This choice reduces equation (6) to 



E(xo,io) 









f ns 


< Hv X 


nvxy 







lis X (lis X dtJ{x^,Q) 



dK 



(14) 



+/3 iiv X J 



lis X at^J(Xs,ts 



dK 



Xs 



where N is the total number of electrons within the interaction area, which, 
presumably, is the same at any time instant. 

Consider a single harmonic component of the time derivative of the source 
current density 

ia;sJ(xs, LUs) exp(ia;sts) = '^'^sJ(xs, <^s) exp ^— i— |xs|^ exp (iujt) , (15) 



where 



t = t. 



1+ pjUs ■ lly) 

1 - /3(no ■ iiv) 
Ixnl 



CCs 



1 + p{n, ■ riy) 
1-/3 



(17) 



Remarkably, frequency relation (16) coincides with the one obtained in [7] for 
the plane-wave case, which did not involve any approximations. We see now 

that even to a purely harmonic but spatially extended source there corresponds 
a whole set of frequencies in the scattered wave, which depend on the particular 
values admitted by iis and in (16). To obtain a more elaborate relation 
between the spatial distribution of the source current density and the spectrum 
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of the scattered wave we take the Fourier transform of the observed electric 
field strength with respect to time. But first we simplify (14) by resorting to 
the effective two-dimensional aperture Ag (see Fig. 1), so that for a purely 
harmonic source (15) we have 



E(nv,to) j {iiv X [iiv X e(ns,a;s)] 



(18) 



+j3 [hv X h(ns,a;s)]}exp (ia;t) dA , 



where e and h are the effective transverse fields with inhomogeneous angular 
distribution. All unimportant coefficients here and below are lumped with the 
constant of proportionality C. Notice that factor exp {iujt) cannot be taken 
outside the integral, since ui depends on the integration variable. The Fourier 
transform of the left-hand-side must be taken with respect to the observation 
time. Whereas, using the following simple manipulation: 



exp I —I 



oo 

(hv, uJo) = j E(nv, to) exp {-iujoto) dU = 

— oo 

oo / I I \ 

-|xo|^ j E fnv,t j exp (-ia;ot) dt. 



(19) 



in the right-hand-side we may take the Fourier transform with respect to 
retarded time t. The observed frequency is denoted ujo- Now we substitute 
(18) in (19) and arrive at 



E 



[riy, uJo) = C j {riv X [riv x e(ns,u;s)] 



(20) 



[nv X h(ns, u^)]} 5{uj - uj^) dA , 



where 5 {uj — uJq) is the Dirac delta-function, which for f3 ^ 0, and ujg ^ 0, can 
be alternatively expressed as 



5 {cu — cuo) — S 



1 + f3{n, ■ Uy) 

c^s^^:^ 
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- (Hs • riv) 



(21) 



Upon substitution of this expression in (20), similarity with the Radon trans- 
form or, more precisely, the Doppler transform [2], [3] becomes obvious. The 
vectorial structure of equation (20) is still quite complex, however. The inten- 
sity of harmonics looks somewhat simpler: 



EE 



(le-Ie)+/3' (ih-Ih) 

- |(nv • Ie)|' - |(nv • Ih)|' + 2/3 (n^ • Re{le x Ih}) , 



(22) 
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where overbar denotes complex conjugation, and 

le = C J e(n„ ujs)6 [cu' - (n^ • n^)] dA , (23) 

nseAs 

h = C I h(n„ a;,)5 [u' - (n, • n,)] dA . (24) 
Splitting these up into components orthogonal and parallel to n^, i.e. 

Ie,h = + lih, (25) 

we finally obtain 

(e-e)= (l?-Ie^)+/3^(lf-Ii^). (26) 

Subsequent transformations involve assumptions about the spatial coherence 
of the source, which are beyond the scope of this paper, but can be found for 
example in [9]. 

Anyway, it is clear that, if we consider two incoherent point sources, then 
the observed field has only two Fourier components corresponding to the fre- 
quencies, which can be easily determined from (16). In other words, we shall 
directly observe an image of our sources in the temporal spectrum of the field. 
Separation between the spectral lines depends on the angular separation of 
the sources and the electron velocity. The higher the velocity, the better the 
attainable spatial (angular) resolution. 



3 Discussion 



Unfortunately, there is no such thing in Nature as a purely monochromatic 
source of type (15). Therefore, even in the preceding simple one-dimensional 
example one should expect broadening of each of the two observed lines pro- 
portional to the broadening of the source's spectral line. In terms of the image 
quality this leads to the loss of spatial resolution - blurring. In addition, spec- 
tral resolution of the dispersing instrument (e.g. a prism or a grating) must 
be taken into account. Imagine two identical dispersing instruments. We place 
one of them on the path between the source and the interaction area, so that 
only one spectral line of width Acus is allowed through. Another dispersing 
instrument is used to investigate the spectral content of the scattered field. 
Then, for the two sources to be resolved, the electron velocity must satisfy 



|Ve| 
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where R = 6i;o/|Aa;o| is, by definition, the resolving power of a spectroscope, 
and A(/9 = | sin — sin (/?2| is the angular separation of the sources (the sources 
and the electron beam are coplanar). 

Another anticipated cause of blurring is the variation of the particle velocity 
inside the interaction area, if the latter could not be made sufficiently small. 
The Fourier transform of (18) with a time- varying velocity will then introduce 
an extra convolution in a;o-domain with the temporal Fourier transform of the 
electron velocity. 

Correspondence between the spectrum and the source for two-dimensional 
distributions is not as direct as in the one- dimensional case. To invert the 
Radon transform spectral data for different directions of the velocity vector 
must be collected, i.e. we have to introduce the mutual rotation between the 
source and the particle accelerator around the axis perpendicular to the plane 
of the source. Subsequently, numerical inversion of the Radon transform must 
be carried out. 

In summary, the proposed rclativistic alternative to aperture-based imaging 
consists of employing a beam of rclativistic electrons as a detector in an optical 
imaging system. A classical treatment of this problem has been presented. A 
more refined analysis would have to take into account: the mutual interaction 
between electrons, the radiation reaction force, and the quantum effects. The 
question of coherence must be investigated in more detail as well. 
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